Superconducting phase transitions in frustrated Josephson-junction arrays on a dice 

lattice 
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Transport measurements are carried out on dice Josephson-junction arrays with the frustration 
index f — 1/3 and 1/2 which possess, within the limit of the XY model, an accidental degeneracy of 
the ground states as a consequence of the formation of zero-energy domain walls. The measurements 
demonstrate that both the systems undergo a phase transition to a superconducting vortex-ordered 
state at considerably high temperatures. The experimental findings are in apparent contradiction 
with the theoretical expectation that frustration effects in the f — 1/3 system are particularly 
strong enough to suppress a vortex-ordering transition down to near zero temperature. The data for 
f — 1/2 are more consistent with theoretical evaluations. The agreement between the experiments 
and the Monte Carlo simulations of a XY model for / = 1/3 suggests that the order-from-disorder 
mechanism for the removal of an accidental degeneracy may still be effective in the f — 1/3 system. 
The transport data also reveal that the dice arrays with zero-energy domain walls experience a much 
slower critical relaxation than other frustrated arrays only with finite-energy walls. 

PACS numbers: 74.81.Fa, 75.10.Hk, 64.60.Cn 



Vortex matters in Josephson-junction arrays (JJA's) 
and superconducting wire networks exposed to a mag- 
netic field have been extensively studied in recent years.~ 
Studies of the systems of a variety of geometries have 
been concentrated especially on the nature of ordering 
transitions and excitations in the ordered state. More 
recently, special attention is given to highly frustrated 
systems which are unable to simultaneously satisfy the 
competing interactions and consequently have the highly 
degenerate ground st at es 1 1 1 1 1 1 2 Understand- 
ing how a highly frustrated system orders or fails to order 
is one of most challenging problems in current condensed 
matter physics. 

The infinite degeneracy of the ground state has been 
suggested in the systems on the triangular, kagome, hon- 
eycomb, or dice lattice with the formation of zero-energy 
domain walls in the ordered vortex state^^^ ' 11 ! 12 The 
common cause to lift the accidental degeneracy is ther- 
mal fluctuations in the continuous phase variable, that 
The difference in free energy of the 
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is, spin waves* 
fluctuations forces many systems to select a particular 
vortex pattern at finite temperatures. This mechanism 
of lifting an accidental degeneracy is often referred to as 
order-from-disorder. 

Frustrated JJA's on the dice (or T3) lattice have at- 
tracted recent attention because of the possibility that 
the order-from-disorder mechanism for the removal of an 
accidental degeneracy may not be effective in the sys- 
tems. Theoretical evaluations^ on the frustrated XY 
model, which most closely represents a JJA in a mag- 
netic field, have demonstrated that for the frustration 
index (defined as the magnetic flux through a rhombus 
tile in units of flux quantum) / = 1/3, frustration effects 
on a dice array are particularly large and the system fails 
to order down to zero temperatures. The stabilization of 
a specific vortex pattern was expected to appear at ex- 
tremely low temperatures only when one goes beyond the 
limits of the XY model by including the magnetic inter- 



actions of currents present in a real JJA. The predicted 
vortex-ordering transition temperature for / = 1/3 is as 
low as 0.01 in units of J /kg, where J is a junction cou- 
pling strength and kg the Boltzmann constant. However, 
it was also expected that the dynamically quenched relax- 
ation of the vortex patterns may prevent an observation 
of the vortex-pattern ordering in experiments. For the 
fully-frustrated dice array with / = 1/2, a commensurate 
vortex pattern was expected to be similarly stabilized at 
T < O.U/fcs by the magnetic interactions of currents in 
a JJA rather than by the fluctuations A 2 - It was argued 
for / = 1/2 that the order-from-disorder mechanism be- 
comes effective only at extremely low temperatures due 
to the reduced contribution to free energy from the fluc- 
tuations as a consequence of a hidden gauge symmetry 
and thus the magnetic interactions of current is a more 
important source for the ordering of vortices. 

The results of theoretical evaluations contradict the 
Monte Carlo (MC) simulations^ of a XY model ex- 
hibiting a periodic vortex pattern for / = 1/3 at T < 
0.2J/ks as a consequence of a phase transition. For 
/ = 1/2, the simulations present a Berezinskii-Kosterlitz- 
Thouless-like (BKT-like) transition to a dynamic glass 
state at T ~ 0.05J/fc s . 

The Bitter decoration experiments 9 on superconduct- 
ing wire networks have also revealed a periodic vortex 
pattern in the / = 1/3 case and no trace of regular vor- 
tex ordering for / = 1/2. For superconducting wire net- 
works, both the magnitude and the phase of supercon- 
ducting order parameter can vary near the temperatures 
of interest. Hence, a quite different situation may de- 
velop in JJA's where the phase fluctuations play a more 
important role. 

This paper presents the results of transport measure- 
ments on dice JJA's with / = 1/3 and 1/2. The results 
show a substantial departure from the theoretical behav- 
ior for / = 1/3. The agreement between the experiments 
and the Monte Carlo simulations for / = 1/3 suggests 
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FIG. 1: (top) Photograph of a dice array used in the ex- 
periment, (bottom) Sample resistance with 30-fiA excitation 
current plotted against frustration at different temperatures 
T = 4.4, 4.25, 4.05, and 3.9 K in descending order. 



that the order-from-disorder mechanism may still be ef- 
fective for ordering in the / = 1/3 system. The data for 
/ = 1/2 are more consistent with theoretical evaluations. 

The experiments were performed on a dice array of 
404x525 Nb/Cu/Nb Josephson junctions. The photo- 
graph in Fig. Q] shows a portion of the JJA sample used 
in this work. Nb islands with a thickness of 0.2 /<m and 
an arm width of 4 /mi were disposed on a 0.3-/mi-thick Cu 
film periodically with the elementary side length of the 
rhombus cell of 16 jxm and a separation between adjacent 
islands of 1.4 /im. The variation of the junction separa- 
tion in the sample was less than 0.1 /im. Measurements 
of the resistance and the TV characteristics were carried 
out by using the standard four-probe technique. The 
sample voltage was measured by a transformer-coupled 
lock-in amplifier with a square-wave current at 23 Hz. 
The single-junction critical current i c at low tempera- 
tures was obtained directly from the IV curve. The i c 
and the junction coupling strength J(= hi c /2e) at high 
temperatures were determined by extrapolating the i c vs 



T data at low temperatures by the use of the de Gennes 
formula^ for a proximity-coupled junction in the dirty 
limit. Additional details of the measurements are de- 
scribed in Ref. 14. 

Figure [1] displays the resistance of the sample with 
30- fi A excitation current plotted against / at four dif- 
ferent temperatures. Local minima are present at low 
order rationals / = 1/6, 1/3, 2/3, and 5/6. A novel 
feature appears at / = 1/2. The anomalous positive 
cusp at / = 1/2 does not disappear even at low temper- 
atures close to the superconducting transition tempera- 
ture T c ~ 3.85 K, determined from IV characteristics 
measurements (to be described below). The f-R curve 
semi-quantitatively reflects the increase in the mean-field 
superconducting transition temperature as a function of 
/. One can thus estimate from the curve the energy cost 
for different vortex configurations. In general, the ground 
state at simple rational / has a simple periodicity^ The 
ground states for other /'s were found to consist of do- 
mains of nearby simple rational /'s i 15 ' 16 ' 17 A finite energy 
required for exciting a domain wall causes a local mini- 
mum of the system energy at the simple rational /'s. If 
the low T c and the formation of a dynamic vortex glass 
proposed for / = 1/3 in theory were the case, a posi- 
tive cusp should appear at f = 1 / 3 . The appearance of a 
distinct minimum at / = 1/3 in the f-R curve indicates, 
against the theoretical prediction^, that the system with 
/ = 1/3 may have a vortex-ordered state with a locally 
minimal energy and a simple periodic structure at consid- 
erably high temeperatures. The positive cusp at / = 1/2 
suggests that the system fails to order within mean-field 
approximations and that excitations of domain walls cost 
an arbitrarily small energy. Apart from the f-R data, the 
T-R data in Fig.[2]of the sample with a smaller excitation 
current of 1 /J.A demonstrate that the f — 1/2 system as 
well as the / = 1/3 system may eventually make a su- 
perconducting transition to an ordered state when the 
temperature is reduced below ^3.9 K. 

The superconductivity for / = 1/2 is verified by the IV 
data plotted in Fig. [31 The distinct activated character of 
the low-temperatures IV curves reveals a superconduct- 
ing state with a long-range phase order as the ground 
state. As displayed more plainly in the d(logF )/d(logJ) 
vs I/I c plot of the same data, the low-temperature ac- 
tivated character changes to the high-temperature resis- 
tive character at T ~ 3.85 K, corresponding to 0.10 in 
units of J/ks- The IV characteristics are quite simi- 
lar to what observed in frustrated square arrays expe- 
riencing a vortex-solid melting transitionj 14 i 17 i 18 i 19 The 
similarity suggests that the superconducting-to-resistive 
transition at T ~ 3.85 K is a melting transition of a 
vortex solid driven by domain walls. Even though the 
f-R data are compatible with the MC simulations^ ex- 
hibiting the dynamic glass state at low temperatures, the 
IV data appear to be more consistent with the theoret- 
ical evaluations^. The positive cusp at / = 1/2 can 
be related with the extreme smallness of the fluctuation- 
induced free energy of zero-energy walls. 
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FIG. 2: Temperature dependence of the sample resistance 
with 1-iiA excitation current for / = 0, 1/3, and 1/2. The 
resistance is plotted in a log scale. 



The IV data for / = 1/3 are shown in Fig. H The 
evolution of the IV data with temperature is similar to 
that for / = 1/2. The IV curves at low temperatures 
have an exponential form, indicative of a vortex-ordered 
state with a long-range superconducting phase coherence. 
With the temperature raised, the low-temperature acti- 
vated behavior crosses over to a resistive form. The IV 
data along with the f-R data indicate that the / = 1/3 
system undergoes a superconducting phase transition 
into an vortex-ordered state with a simple periodicity 
at T ~ 4.15 K (= 0.2J/fc B ). The T c ~ Q.2J/k B is 
comparable to those of other frustrated arrays without 
any accidental degenerac y 14 i 17 i 19 The occurrence of a 
superconducting phase transition into a vortex-ordered 
state at such a high temperature is in apparent contra- 
diction with the theoretical argument 1 ^ that frustration 
effects for / = 1/3 are particularly large enough to pre- 
vent the system from developing any order down to near 
zero temperature. However, it is in agreement with the 
MC simulations 10 of a XY model and compatible also 
with the decoration experiments 9 - on wire networks. The 



FIG. 3: Evolution of IV curves with temperature for / = 1/2. 
The 7 vs V isotherms differ by a temperature interval of 0.1 
K. The panel at the bottom shows the slope of the isotherms 
as functions of I/I c . The dashed line is drawn to show where 
the superconducting phase transition occurs. 

agreement between the experiments and the simulations 
may imply that the removal of an accidental degeneracy 
for f — 1/3 can be explained within the frame of a XY 
model or that the order-from-disorder mechanism may 
still be effective in the / = 1/3 system. 

Another interesting feature of our data is that the dy- 
namic critical exponent z of the systems determined from 
the IV data is unusually high. For a continuous su- 
perconducting transition, the IV data satisfy a simple 
power-law IV relation, V ~ I z+1 , at T c . 20 One can thus 
obtain z directly from the straight log/ vs logV^ isotherm 
at T c . z of the dice arrays determined from the IV data 
is ~ 1.7 for both / = 1/3 and 1/2, which is much larger 
than those (z < 1) of other frustrated array s 14 ' 17 ' 19 ex- 
periencing a vortex-lattice melting transition driven by 
finite-energy domain walls. The exponential IV charac- 
teristics at T < T c do not allow one to relate the large 
z to the possible occurrence of a BKT-like transition. 
The large z demonstrates that the critical relaxation in 
the dice arrays with an accidental degeneracy is much 
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slower than in other frustrated arrays without such an 
accidental degeneracy. It may imply that the motion of 
zero-energy domain walls is not so effective in relaxing 
the system near T c as that of finite-energy walls is. 



FIG. 4: Evolution of the IV curves with temperature for 
/ = 1/3. 



To summarize, transport measurements were carried 
out on dice JJA's with / = 1/3 and 1/2 which possess, 
within the limit of the XY model, an accidental degen- 
eracy of the ground states as a consequence of the for- 
mation of zero-energy domain walls. The measurements 
demonstrate that both the systems undergo a phase tran- 
sition to a superconducting vortex-ordered state at con- 
siderably high temperatures. The experimental findings 
are in apparent contradiction with the theoretical expec- 
tation that frustration effects in the f — 1/3 system are 
particularly strong enough to suppress a vortex-ordering 
transition down to near zero temperature. The data 
for / = 1/2 are more consistent with theoretical eval- 
uations. The agreement between the experiments and 
the MC simulations of a XY model for f — 1/3 suggests 
that the order-from-disorder mechanism for the removal 
of an accidental degeneracy may still be effective in the 
/ = 1/3 system. The transport data also reveal that 
the dice arrays with zero-energy domain walls experience 
a much slower critical relaxation than other frustrated 
arrays only with finite-energy walls. 
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